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Abstract 

This report presents computations for the Type !V shock-shock 
interference flow under laminar and turbulent conditions using unstructured 
grids. Mesh adaptation was accomplished by remeshing, refinement and mesh 
movement. Two two-equation turbulence models were used to analyze 
turbulent flows. The mean flow governing equations and the turbulence 
governing equations are solved in a coupled manner. The solution algorithm 
and the details pertaining to its implementation on unstructured grids are 
described. Computations were performed at two different freestream Reynolds 
numbers at a freestream Mach number of 11. Effects of the variation in the 
impinging shock location are studied. The comparison of the results in terms of 
wall heat flux and wall pressure distributions is presented. 

Nomenclature 


A Jacobian matrix 

a speed of sound, m/sec 

B matrix associated with the viscous fluxes 

BT boundary term 

C constant 

Cf local skin friction coefficient 

C p specific heat at constant pressure, J/kg-K 

C v specific heat at constant volume, J/kg-K 

E total energy, J/kg 

F inviscid flux 

f damping function 

G viscous flux 

h specific enthalpy, J/kg 


k turbulent kinetic energy, J/kg 

M Mach number; mass matrix 

N element shape functions; number of sides 

n direction cosines 

P production terms 

p pressure, N/m 2 

Pr Prandtl number 

q heat transfer rate, MW/m 2 

q 0 stagnation heat flux in the undisturbed flow, MW/m 2 

r radius of the cylinder, in 

R eigenvector matrix 

Re Reynolds number 

S source terms 

T temperature, K 

t time, sec 

U conservation variable 

u component of velocity in x-direction; velocity component 

in indicial notation, m/sec 
v component of velocity in y-direction, m/sec 

x,y coordinates of the reference frame 

GREEK SYMBOLS 
a constant 

5 Kronecker delta; length of an element side; 
boundary layer thickness 

e turbulent dissipation rate, m 2 /sec 3 

6 eigenvalue limiter 

y ratio of specific heats 

r edge of an element 

A eigenvalue 

p viscosity, N-s/m 2 

v kinematic viscosity, m 2 /sec 

0 angle 

p density, kg/m 3 

<j normal stress, N/m 2 

t shear stress, N/m 2 

co specific dissipation rate, 1/sec 

£2 area of an element, m 2 

SUBSCRIPTS 

i,j,k indices of notation 

inf freestream conditions 

e element 

L left value; lumped value 
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n normal direction 

o stagnation value 

peak maximum 
s side 

T turbulence quantity 

w wall 

x,y coordinate axes 

oo freestream 

SUPERSCRIPTS 

* current value 

m time step index 

+ near wall non-dimensionalization 

Introduction 

Intense aerothermal loads that result from shock-shock interferences 
(especially Type III and Type IV interferences) on a hypersonic vehicle are 
critical for design purposes. Computational methods are being developed at 
the NASA Langley Research Center to accurately predict these loads. These 
efforts have led to the development of LARCNESS (acronym for Langley 
Adaptive Remeshing Qode NaviEr Stokes Solver) which is a software package 
that includes a flow solver based on a point-implicit, upwind algorithm, and an 
unstructured mesh generator based on an adaptive remeshing technique. 
LARCNESS has been used successfully in resolving complicated flow features 
that result from shock-shock interferences under laminar flow conditions. The 
details of the algorithm and the computational results are found in Reference 1. 

In a typical Type IV shock-shock interference problem a supersonic jet is 
produced which impinges on the body as shown in Fig. 1. This impingement 
causes highly localized peaks of wall pressures and wall heat fluxes. 
LARCNESS was used to analyze this problem computationally using 
unstructured grids {Ref. 1}. There is some experimental evidence to show that 
the shear layers surrounding the supersonic jet are turbulent during 
experiments {Ref. 2}. This indicated a need for implementation of a turbulence 
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model in LARCNESS to investigate the turbulence effects, and this effort 
addressed this problem. 

Algebraic turbulence models (mixing-length type) have been 
successfully used in compressible turbulent flow computations (for example, 
see Ref. 3). However, their application is limited to wall bounded flows using 
structured grids. They also lack generality in the context of application to 
complicated geometries and complex flow fields such as shock waves, jets and 
shear layers. On the other hand, the two-equation models are not limited to wall 
bounded flows, and offer the flexibility of dealing with complex flow features 
using unstructured grids. In general, these models introduce two turbulence 
parameters, one for the turbulence kinetic energy, k, and the other for the 
turbulence length scale determining parameter (for example, the turbulent 
dissipation rate, e ). These turbulence parameters which define the turbulent 
viscosity in the flow field are governed by two partial differential equations that 
must be solved together with the mean flow equations. Although there are 
several higher-order turbulence models like the Reynolds-stress models, they 
are still in the developmental stage, leaving the most popular two-equation 
models for practical applications. 

The two-equation model was first developed by Jones and Launder {Ref. 
4}. Several modifications to this model have been proposed recently, and 
concise reviews on these models may be found in References 5, 6 & 7. Among 
the many two-equation models available, only a few are capable of modeling 
high speed compressible turbulence. One such model is the compressible k-e 
model with near-wall damping functions developed by Zhang, et al. {Ref. 8}. 
This is the most widely used two-equation model in the literature, and it has 
been very successful for a large variety of different flow situations, including free 
shear layers. However, the k-e model suffers from several shortcomings. 
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Some of these shortcomings include lack of sensitivity to adverse pressure 
gradients, problems related to numerical stiffness in the near-wall regions and 
in some cases underprediction of wall heat fluxes. {Ref. 9 & 10}. 

On the other hand, another very successful two-equation turbulence 
model, the k-co model originally developed by Wilcox {Ref. 11} is robust and 
accurate in predicting near wall turbulence even under adverse pressure 
gradients. The major shortcoming of the k-co model is its high sensitivity to 
freestream conditions and hence its inability to predict free shear layers 
accurately. A combination of the k-e model and the k-co model appears to 
eliminate each other's shortcomings and give rise to a better turbulence model 
{Ref. 12}. Recently several zonal approaches combining these models were 
successfully attempted in solving high speed flow problems {Ref. 12 & 13} and a 
similar approach is used in the current work. Since, in the present code, 
structured grids are used near the wall and unstructured grids elsewhere, the k- 
co model is used in the structured grid portion of the mesh and the k-e model is 
used in the unstructured grid portion of the mesh. 

Some work has been reported in the literature regarding the 
implementation of two-equation turbulence models for compressible flow in the 
context of unstructured grids. Marcum and Agarwal {Ref. 14}, implemented the 
k-e model in a Galerkin finite-element procedure for solving the Reynolds- 
averaged mean flow equations. Mavriplis {Ref. 15} implemented a 
compressible version of the k-e model in a multi-stage Runge-Kutta time- 
stepping scheme using a multigrid strategy on unstructured grids. In both of 
these papers however, the mean flow equations and the turbulence equations 
were solved in an uncoupled manner. Morrison {Ref. 16} implemented the k-e 
model proposed by Zhang, et al. {Ref. 8} in a cell-centered finite-volume 
technique, using the Roe's flux difference splitting technique in a coupled 
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manner (solving all the governing equations simultaneously) on structured 
grids. In the present work, unstructured grids are employed, and both methods 
of solving (coupled and uncoupled methods) the governing equations are 
investigated. 

The purpose of this report is to present the laminar and turbulent flow 
computations using two-equation turbulence models in association with an 
unstructured grid solver for the Type IV shock-shock interference problem. First 
the turbulence models, along with some details regarding its implementation 
into the LARCNESS code, are presented. The computational results for the 
shock-shock interference at two different freestream Reynolds numbers in a 
Mach 1 1 flow over a 3 inch cylinder are presented. A comparison of the laminar 
flow computations with the turbulent flow computations are made. Conclusions 
are drawn based on these studies. 


Governing Equations and the Turbulence Models 

The basic governing equations for the mean flow are the Favre-averaged 
Reynolds equations along with Boussinesq approximation and are given as 
follows: 
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("tilde" over-bar represents Favre mass averaged quantity, and plain over-bar 
represents the time averaged quantity). In the above equations, jtr is the 

turbulent viscosity, Pr 7 is the turbulent Prandtl number and k is the kinetic 
energy of turbulence given by 


k=ptf'iu"j/2p (4) 

The k-e model 

In the k-e model, the turbulent viscosity is expressed in terms of turbulent 
kinetic energy k, and the turbulence dissipation rate e as follows: 


— C^pk 2 
fi T = -~ 


( 5 ) 


where Cp is a constant and the dissipation rate e is expressed as follows: 
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This model introduces two more unknowns in the original set of governing 
equations, i.e„ the turbulent kinetic energy (k) and the turbulence dissipation 
rate (e), which are governed by the following partial differential equations 


expressed in the conservation form : 
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In the above equations, Pk represents a production term given as follows 


Pk=- pL/'lLf)^ = H T 
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where CiC 2 , c>k, o e , and a are constants, and M r is the local turbulent Mach 

number defined in terms of turbulence kinetic energy and local speed of sound 
a, as follows: 



(10) 

The constants used in Eqn. (7) and (8) are given below 

(11) 

C 1 = 1.44 


C 2 = 1.92 


(77= 10 


O e = 1.0 


a = 0.5 



The k-co mode) 

In the k-co model, the turbulent viscosity is expressed in terms of turbulent 

kinetic energy k, and the specific dissipation rate co as follows: 

— pk 

VT=-bz- ( 12 ) 

CO 

where co is defined as 

= (13) 

Cfx k 

The governing equation for the turbulent kinetic energy k in this model is the 
same as the governing equation in the k-e model except the dissipation terms 
are expressed in terms of the specific dissipation rate co, and given as follows: 
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In order to incorporate the low turbulent Reynolds number effects (that are 

present in the laminar sublayer region near a solid wall) into the turbulence 

model, a set of damping functions as suggested in Ref. 11 are implemented in 
the current model. 

The cross-diffusion term 

Since the turbulence dissipation rate, i is related to the specific 
dissipation rate 5, equation (8) can be written in terms of £ as follows- 


(17) 




+ 2 aw5/£ + _£_j 
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Comparing Eqns . (8) and {17) it can be easj|y notjced ^ ^ djfference 

between the e equation expressed in terms of 5 and the 5 equation is the cross- 
es, on term that appears at the end of right-hand-side of Eqn. (17). In other 
words dropping of the cross-diffusion term from the transformed i equation 
results exactly in the governing equation for S. Since the governing equation 
for At ,s the same in both the models many of the shortcomings of the ** mode, 
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are attributed {Refs. 12 & 17} to the presence of this cross-diffusion term in the e 
governing equation . Conversely, the shortcomings of the k-co model are 
generally attributed to the absence of this cross-diffusion term in the 5 
governing equation. In the zonal approach, the formulation is made according 
to the k-co model and the cross-diffusion term is introduced in the co governing 
equation in the regions where the k-e model is employed through a blending 
function. In the current formulation, co is used as a length scale turbulence 
parameter and the cross-diffusion term is introduced into the co governing 
equation to recover the k-e model wherever needed. 


Implementation in LARCNESS Code 

The set of governing equations are solved using an upwind, point- 
implicit, cell-centered scheme based on Roe's averaging. The point-implicit 
time marching procedure for the Navier-Stokes equations is given in Ref. 1. 
Extension of this solution procedure for Reynolds averaged mean flow 
equations was accomplished by including additional terms in the mean flow 
equations. In the current approach the four mean flow equations and the two 
turbulence governing equations are solved in a coupled manner, as proposed 
in Ref. 16. 


The details of the solution procedure for the coupled set of governing 
equations are as follows. For the sake of simplicity, over-bars are eliminated in 
the following discussion. The set of governing equations can be written in the 
conservation form as 


iM + ^ = .aM +{s) 


( 18 ) 


d t d X, d Xj 

where {U} is the vector of conservation variables, {F} is the inviscid flux vector, 
{G} is the viscous flux vector, and {S} is the vector containing the source terms 
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(source terms exist only in the turbulence equations comprising of production 
and dissipation or destruction terms). These vectors are given as follows: 
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where a xx , a yy and r xy are the stress components, D kx , D ky , D mx , and D ojy 
are the diffusion terms in the turbulence equations, and q x and q y are the 
components of the heat-flux vector. A calorically perfect gas is assumed. 

The stress components are given by 


°*x = § (fl + flj) 


^ 3 U 3 V 


d X 3 y 


+ 3 pk 


T xy- ?yx- {H + /It) 


3 U 3 V 
3 y 3 x 


and 


°yy = f(M + ^r) 


2 3 v 3 u 


I 3 y dx 

The components of the heat-flux vector are given by 




and 
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The diffusion terms in the turbulence equations are given as follows: 
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The solution domain is discretized spatially into an unstructured grid of 
triangular and quadrilateral elements, and the set of governing equations is 
solved subject to appropriate initial and boundary conditions. The governing 

equations for each element can be expressed in an integral manner as given 
below. 
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where D e is the area of the element. Assuming a piecewise constant 
distribution of the unknowns U and linearizing l/with respect to t, Eqn. (29) can 
be rewritten as shown below. 


i Ue^-jUs)" 1 _ 1 
At Q e 



a{ G)j 3(F) A 
l dx i ' ' ix i ] 


d£2+ {S)dQ 


(30) 


where U e represents element level unknowns and m denotes the time level. 

By application of the divergence theorem to the flux terms on the left-hand-side, 
Eqn. (30) can be written as follows: 

l ««■«' <3» 

</ J. Q 

where r e is the element boundary and the subscript n denotes the flux normal 
to the element sides. 

Assuming the viscous fluxes to be a constant along the sides of an 
element and the source terms to be constant within the element, Eqn. (31) can 
be written as 

{AU e ) =M. £ 5 s [{G n r 1 -{F„}' m1 ] s +r2 e (S e }' m - 1 (32) 

&e \ S=1 I 

where S s is the length of the side under consideration, and nside denotes the 
total number of sides in the element. In Eqn. (32) G n and F n represent the 
viscous and inviscid fluxes respectively, and S e denotes the source terms. 
Computation of these quantities is discussed in the following section. 

Treatment of the Fluxes 

The inviscid flux vector F n appearing in Eqn. (32) is replaced by a 
numerical flux given by 



1 3 



( Fn)=l[{FR) + (F L ) - [|4] ({U R )-{U L })] (33) 

where subscripts R and L denote the right and left values with respect to the 
side s, and the term [\A\] {{Ur} -{Ul}) is the first-order dissipation term as 
suggested by Gnoffo {Ref. 18}. The matrix [A] is the Jacobian of the inviscid flux 
vector evaluated at the Roe's average state {Ref. 19}, and has the property 

{Fr)-{Fl)=[A]({U r )~ {U L }) (34) 

This matrix can be factored as 

[4l = [fl]" 1 U] [H I (35) 

where [A] is the diagonal eigenvalue matrix of [A] and [R] is the corresponding 

eigenvector matrix of [>4]. The complete form of the Jacobian matrix [>4] and the 
matrices [R] and [R]~ 1 are given in Ref. 16. 

The matrix [ \A\ ] is defined as 

[H] = [fl]' 1 [ W][fl] (36) 

where [ IAI ] is the diagonal absolute eigenvalue matrix of [A]. Matrix [>4] has 

six eigenvalues A/, and the minimum value for the eigenvalues A is restricted 
according to Harten {Ref. 20} such that 
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0.5 







<eu| 


( 37 ) 


where ex is the eigenvalue limiter. 

The determination of the numerical viscous fluxes involves evaluation of 
the first derivatives of the components of velocity, the temperature, the 
turbulence kinetic energy, and the turbulence dissipation rate at the nodes. In 


the cell-centered scheme discussed above, the primitive variables are assumed 
to be piecewise constant. That is, they are assumed to be constant for an 
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element. In order to compute the first derivatives at the nodes, a bilinear 
variation of these primitive variables is assumed for quadrilateral elements, and 


a linear variation is assumed for triangular elements. Under these assumptions, 
the first derivative of a variable, T, for example, at any point in an element can 
be written as 


37 

dx 


= [N] 


where [N] are the shape functions and 



(38) 

is the nodal first derivative vector 


of the element. By multiplying left-hand-side of Eqn. (38) by {N} and on 
integration by parts we obtain 



where n are the direction cosines. By multiplying the right-hand-side of Eqn. 
(38) by {N} and on integration we obtain 



{A/} [A/] 




{N}[N]dQ 



(40) 


= M 

where [MJ is the lumped mass matrix obtained by diagonalizing the consistent 
mass matrix or the capacitance matrix. By substituting Eqn. (40) in Eqn. (39) 
and rearranging, we obtain the following expression: 




nT{N]dr - 






(41) 
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To compute the first derivative at a particular node /, an assembly operation is 
performed where contributions from all the elements (ne) that have /as a 
common node are summed up. This operation results in the following relation 
for the first derivative at node /. 


Bx j m 


BT xi 



(42) 


where is the coefficient of the global (assembled) lumped mass matrix for 
node /, BT X j is the assembled boundary term (needs to be evaluated only at the 

domain boundary) surrounding node i, and T e is the temperature of the element 
e. 


The source terms in the turbulence equations are computed at the 
element level. The velocity derivatives needed at the element level to compute 
these source terms are obtained by averaging the nodal derivatives of the 
element. The cross-diffusion terms required for the k-e model are also 
computed at the element level, similar to the source terms. 

Point-Implicit Time Stepping Technique 

Substituting Eqn. (33) for the inviscid numerical flux in Eqn. (32) and 
rearranging, we have 

i nside 

{Au e } =m [ .±2 5 s [(F fl r i +{F L r i -[|/iD'™- i ({t; n r i -{^} mfi )] s 

&e d s= 1 

nside ( 43 ) 

+2 ^(Gsr'-kG^r 1 } 

S=1 

where {G s } represents the numerical viscous flux, vector for the side s and the 
subscript e represents the element under consideration. The time integration 
scheme as shown above is a fully implicit technique and leads to solving the 
entire global system of equations. Equation (43) can be linearized by 
employing a Gauss-Seidel point relaxation scheme which results in a point- 
implicit scheme where the set of governing equations is solved at each element 



as given in {Ref. 1}. The point-implicit scheme is represented in the following 
equation. 
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Q, e 


aside 

<4X 
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8, [{ F r )\{ F t } m M({ U R )'-{ I4)i s 


aside 

+ 1 Si{G s ) m ] + Q e {S e n (44) 

s= 1 

where superscript * refers to the latest available value. Matrix [B s ] represents 

the implicit contribution from the viscous fluxes and the matrix [— ] denotes the 

dU 

Jacobian of the source terms in the governing equations for the turbulence 
variables. 


LARCNESS Computations 

This section describes the application of LARCNESS to the problems of 
shock-shock interferences in hypersonic flow. Comparison of laminar flow 
computations with that of turbulent flow are made wherever applicable. The 
turbulence models used in this work do not predict transition and in all the 
applications discussed in this section, no account was made for the transition 
from laminar to turbulent flow. For the application problems, an initial solution is 
obtained without the source terms in the turbulence governing equations and 
these terms are introduced later into the equations to obtain a final converged 
solution. The freestream turbulent viscosity is always assumed to be of the 
same order of magnitude as the freestream laminar viscosity and the freestream 
turbulent kinetic energy is assumed to be 0.1% of the freestream mean flow 
kinetic energy. The validation of the code was performed by two test cases and 
the results are presented in {Ref. 9}. 

The computer program for the solution algorithm is fully vectorized on the 
CRAY computers and requires a memory of 11.8 MW for the turbulent flow 
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computations on a 200,000 element mesh compared to 9.4 MW of memory for 
laminar computations on the same mesh. On a CRAY-C90 computer, 
LARCNESS code took 8.44E-06 CPU seconds/element/iteration for laminar 
flow computations (for a 35,000 element mesh) compared to 1.07E-05 CPU 
seconds/element/iteration for the turbulent computations on the same mesh. It 
is interesting to note that the increase in the CPU time for the turbulent 
calculations is only 20% more than that for the laminar calculations. 

Type IV Shock-Shock Interference in a Mach 11 Flow 

The wall pressure and heat flux amplification is a maximum for the Type 
IV shock-shock interference relative to the rest of the interferences as 
mentioned in Ref. 21. A schematic of the Type IV interference is shown in Fig. 1. 
In a typical Type IV interference, the impinging shock interferes with the near 
normal portion of the bow shock resulting in a supersonic jet that impinges on 
the body. This supersonic jet terminates through a normal shock before it 
reaches the body, creating a small stagnation region where the heat fluxes are 
the maximum. The supersonic jet is surrounded by shear layers, the nature of 
which (either laminar or turbulent) plays a role in augmenting the wall heat 
fluxes {Ref. 22}. To investigate the effects of turbulence this problem was run 
with both the laminar and turbulent options under two sets of freestream 
conditions as given in Table A. These freestream conditions are for air which is 
considered as a perfect gas with ratio of specific heats, y, of 1 .4. The first set of 
conditions, referred to as condition A (throughout this report), is at a freestream 
Mach number of 11.36 and freestream Reynolds number of 6500 based on the 
cylinder radius, and is assumed to correspond to laminar shear layer flow 
conditions. The second set, referred to as condition B (throughout this report), is 
at a freestream Mach number of 11.03 and freestream Reynolds number of 



414,000 and is assumed to correspond to turbulent shear layer flow conditions. 
These freestream condition sets correspond to experimental test conditions to 
be reported at a later date by R.J. Nowak***- For both the sets of freestream 
conditions, the wall temperature was fixed at 530°R. Since the freestream 
conditions played a major role in the computations, results from LARCNESS 
computations are presented in two separate subsections for the two condition 
sets in the following discussion, 

Freestream condition set A 

This set of freestream conditions corresponds to a lower Reynolds 
number as mentioned above. Before the shock-shock interference problem 
was attempted, an undisturbed flow (without the shock-shock interference) 
under the same conditions over a 1.5 inch radius cylinder as shown in Fig. 2 
(problem statement) was analyzed using the LARCNESS code. The theoretical 
stagnation heat flux, q 0 , was computed from the Fay and Riddel theory {Ref. 23} 
for the undisturbed flow without the shock-shock interference . The wall heat 
flux distribution q from the LARCNESS computations non-dimensionalized with 
respect to the theoretical stagnation heat flux q 0 is shown in Fig. 3. The 
computational stagnation heat flux was within 2% of the theoretical stagnation 
heat flux. 

For the shock-shock interference flow problem as shown in Fig. 4, 
successive mesh adaptation was performed until the last two meshes 
essentially gave the same result. These meshes consisted of: (a) highly 
stretched structured grids next to the wall, and (b) solution adapted unstructured 
grids everywhere else. The unstructured portion of the mesh adaptation was 
accomplished by a combination of three different adaptive techniques. The first 


Dr. R.J. Nowak can be reached at NASA LaRC by phone at (804) 864-1341 or by mail at M/S 
408A 
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technique is the adaptive remeshing where a new mesh is created based on 
the solution from a previous mesh {Ref. 24}. The adaptive remeshing did not 
produce the level of refinement that was required in the shear layers. Therefore 
the second technique of subdividing the elements in the shear layer region was 
employed in association with the remeshing technique. To further enhance the 
quality of the mesh, a third technique known as adaptive mesh movement was 
used where the nodes in the mesh are moved closer to each other wherever 
higher levels of mesh refinement was desired. This combination of mesh 
refinement techniques applied to high speed flows is described in {Refs. 9 & 
25}. The final mesh consisted of 210,195 elements and 119,959 nodes. The 
size of the smallest element was 4X1 O' 7 inches. A portion of this final mesh 
near the stagnation region is shown in Fig. 5. Convergence criterion was set at 
3 or 4 orders of magnitude reduction in all the residuals corresponding to the 
mean flow variables and less than 1% change in the peak heat flux value over a 
specified period of time {at least 1000 iterations). The convergence in terms of 
peak heat flux history for the last 3000 time iterations on the final mesh for 
condition set A under laminar flow conditions is shows in Fig. 6. Table B 
includes the coordinates of the two triple points and the transmitted shock 
lengths for both the laminar and turbulent cases corresponding to freestream 
condition set A. 

The total velocity contours, pressure contours, and first derivative of 
density contours near the stagnation region, from the LARCNESS computations 
on the final mesh, are shown in Fig. 7, 8, and 9. The x and y axes shown in 
these figures are reference axes parallel to the coordinate axes X and Y (the 
origin of these reference axes is not shown at the center of the cylinder which is 
the origin for the coordinate axes). These contours show the intricate details of 
the flow field the shear layers, the supersonic jet and the jet terminating shock. 
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The jet terminating shock location as shown in these contour plots suggests that 
this case does not represent the peak Type IV case where the jet terminating 
shock is parallel to the cylinder wall and hence the supersonic jet impingement 
is normal to the wall. Turbulent flow computations were performed using the 
LARCNESS code and the total velocity contours near the stagnation region 
from both the laminar and the turbulent computation are compared in Fig. 10. 
The shear layers in the turbulent case appear to be thicker than that of the 
laminar case. The wall heat flux distribution from the turbulent computations 
was compared with that of laminar computations and is shown in Fig 11 A. 
Although there was a shift in the location of the peak heat flux (about 2°), there 
was no significant difference between the magnitudes of the peak heat flux in 
both the cases. The reason that the wall heat flux distributions are not very 
different from each other is that the freestream Reynolds number is very small 
and the influence of the turbulence effects on the flow field are negligible. 
Similarly the wall pressure distributions from the laminar and the turbulent 
computations is shown in Fig 11b. A trend similar to that of the wall heat flux 
distributions was observed for the wall pressure distributions. To investigate the 
influence of the impinging shock location in the inlet plane of the computational 
domain, the original location of the impinging shock was lowered by 0.055 
inches and laminar computations were performed for this new impinging shock 
location. The comparisons of the wall heat flux distributions and the wall 
pressure distributions from the laminar computations for the two impinging 
shock locations are shown in Fig. 12. For this freestream condition set, 
although there was a difference in the locations of the peaks in terms of wall 
pressures and wall heat fluxes as expected, there were no major differences in 
the magnitudes of these two quantities for the two impinging shock locations. 
Table C contains the location of the peaks and values of these peaks (in terms 



of degrees) of wall pressures and wall heat fluxes for both laminar and turbulent 
cases for this freestream condition set A. The magnitudes of the peak pressure 
in Table C are non-dimensionalized with respect to p Q , the theoretical 
stagnation pressure of 0.5159 psia. The magnitudes of the peak heat flux in the 
table are non-dimensionalized with respect to q 0 , the theoretical stagnation heat 
flux of 20.26 Btu/ft 2 -sec. 

Freestream condition set B 

This freestream condition set corresponds to a higher Reynolds number 
and is assumed to be under turbulent shear layer conditions. To understand 
the effects of turbulence on the flow field, this case was run under laminar as 
well as turbulent conditions. The most refined mesh for the turbulent case 
consisted of 245,873 elements and 156,432 nodes. The size of the smallest 
element next to the wall in the normal direction is 4x1 O' 8 inches. A similar level 
of refinement was obtained on the mesh for the laminar case. The convergence 
criterion was the same as that mentioned in the discussion under freestream 
condition set A. The contours of pressure, and the first derivative of density are 
shown in Figs. 13, 14. Again in these contour plots, the axes x and y represent 
the reference axes only and their origin is not at the center of the cylinder. 
Figure 15a show the total velocity contours in the entire flow field. Figure 15b 
shows the total velocity contours near the stagnation region for the laminar case 
and Fig. 15c shows the total velocity contours for the turbulent case. One may 
observe from this figures 15b and 15c that the turbulent shear layers are thicker 
than the laminar shear layers. To obtain an estimate of the shear layer growth 
and thickness, a normal cut parallel to the y-axis was made through the shear 
layer in the region between the 1st and the 2nd triple points (see Figure 1 and 
4). The shear layer thickness at a given x-location is defined as the y-distance 
between the upper and lower points in the shear layer where u-velocities are 
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97.5% of the upper and lower local freestream velocities respectively. At x = 
-2.076 inches from the center of the cylinder, the shear layer thickness for the 
turbulent case was 0.3134 inch; and the shear layer thickness for the laminar 
case was 0.2663 inch (about 15% thinner). However, the shock-shock structure 
including the triple points was farther from the body for the turbulent case 
compared to the laminar case. In order to compare laminar and turbulent shear 
layer thickness at the same distance from the respective 1st triple points, the 
laminar shear layer thickness was calculated at x = -1.973 inches, which is the 
distance from the 1st triple point for the above turbulent shear layer thickness. 
The shear layer thickness for the laminar case at this location was 0.2015 
inches which is significantly thinner (about 67%) than the turbulent case. 
Similar to the case of freestream condition set A, the jet terminating shock is not 
parallel to the wall indicating that this case does not represent the peak Type IV 
location. Therefore, the shear layer effects on the stagnation region are 
predominant. Fig. 16 shows the streamline contours superimposed on the total 
velocity contours for the turbulent case. These contours show that a major 
portion of the flow going towards the stagnation point is coming from the 
supersonic jet region and the upper shear layer. 

The wall heat flux distributions from the laminar and turbulent 
computations are shown in Fig. 17a. The peak heat flux distribution from the 
turbulent computations is significantly higher than peak heat flux from the 
laminar computations. The location of the peak in the turbulent computations is 
also different from the laminar computations and this difference is again 
attributed to the shift in the shock structure as explained above. The 
comparison of the wall pressure distributions for the laminar and the turbulent 
cases is shown in Fig. 17b. The peak pressure in the turbulent case seems to 
be lower than that of the laminar case. Two different sets of computations were 
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performed for two different impinging shock locations. Fig. 18 shows the 
comparison of the wall heat flux distributions and the wall pressure distributions 
for the two impinging shock locations and there is a clear evidence of high 
sensitivity of the peak pressure and peak heat flux to the location of the 
incoming impinging shock. Table E contains the location of the peaks (in terms 
of degrees) of wall pressures and wall heat fluxes for both the cases for this 
freestream condition set B. The magnitudes of the peak pressure in Table E are 
non-dimensionalized with respect to p 0 , the theoretical stagnation pressure of 
7.005 psia. The magnitudes of the peak heat flux in the table are non- 
dimensionalized with respect to q 0 , the theoretical stagnation heat flux of 19.78 
Btu/ft 2 -sec. 


Concluding Remarks 

A two-equation compressible turbulence model is implemented in a 
point-implicit, cell-centered , upwind algorithm applicable to unstructured grids. 
The current turbulence model based on a zonal approach involves k, the 
turbulent kinetic energy and co, the specific turbulence dissipation rate as the 
two turbulence variables. This model results in the standard k-co model in the 
boundary layer next to the wall and results in the standard k-e model away from 
the wall by inclusion of a cross-diffusion term in the governing equation for the 
specific dissipation rate. The governing equations for the turbulence variables 
are solved simultaneously with the equations governing the mean flow. For the 
problem of a 2-D hypersonic flow (Mach 11) over a circular cylinder, the peak 
heat flux from the laminar flow computations agreed well with the peak heat flux 
computed from the Fay and Riddel theory. 

For the Type IV shock-shock interference in a hypersonic flow over a 
cylinder, two sets of freestream conditions were chosen with two different 
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freestream Reynolds numbers. The flow conditions at the higher Reynolds 
number freestream conditions are assumed to be turbulent. For the first set of 
freestream conditions, the results from the laminar and turbulent computations 
did not differ significantly from each other. Highly refined solution adaptive 
meshes were used in these computations to capture the intricate flow details 
such as the shear layers, the supersonic jet, and the jet terminating shock. The 
sensitivity of the augmented pressures and heat fluxes due to the variation in 
the location of the incoming impinging shock seem to be less significant for this 
freestream condition set. 

For the second set of freestream conditions with a higher freestream 
Reynolds number, there were significant differences in the results from the 
laminar and turbulent conditions. The velocity contours showed thicker shear 
layers in the turbulent computations relative to that of laminar computations. 
The wall heat flux distribution from the turbulent computations showed a 40% 
increase in the value of the peak heat flux compared to that of laminar 
computations. The peak wall pressure in the turbulent computations was lower 
relative to the peak in the laminar computations. 

The current work is one of the first attempts to apply two-equation 
turbulence models to the shock-shock interference problem. As the freestream 
Reynolds number goes higher, the turbulence effects seem to influence the flow 
field stronger and result in higher wall peak heat fluxes and wall peak 
pressures. 
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Table A. Freestream conditions used for the Type IV interference analysis 


Variable 

freestream condition set A 

freestream condition set B 

Velocity 

8366.87 ft/sec 

5018.3 ft/sec 

Mach no. 

11.36 

11.03 

Reynolds no. 

6636 

413,974 

Static 

Temperature 

225. 8°R 

83.20° R 

Static 

Pressure 

3.051 E-03 psia 

4.396E-02 psia 


Freestream Reynolds number is based on the cylinder radius of 38.1 mm (1 .5 in) 

*** Flow conditions below the oblique shock were determined based on oblique shock 
relations with a flow deflection angle of 10° and a shock angle of 14.05° 
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Table B. Coordinates of the triple points and transmitted shock lengths for freestream 
condition set A for the impinging shock location 2 


Case 

1st triple point 
coords, (non-dim) 

2nd triple point 
coords, (non-dim) 

Transmitted shock 
length, in. 

1) Laminar 

x = -1 .72069 
y = -0.23732 

x = -1.10384 
y = -0.53386 

1 .02664 

2) Turbulent 

x = -1.70228 
y = -0.22691 

x = -1. 11194 
y = -0.51852 

0.98765 


*** 


All the coordinates are non-dimensionalized with respect to the cylinder radius of 
38.10 mm (1.5 in) 




Table C. Locations and magnitudes of the pressure and heat flux peaks for freestream 
condition set A 


Case 

Location of the 
pr. peak (deg) 

Magnitude of the 
pr. peak 

Location of the 
h.f. peak (deg) 

Magnitude of the 
h.f. peak 

1) Laminar 
(Impinging 
shock loc.1) 

-32.00 

17.00 

•31.63 

23.22 

2) Laminar 
(Impinging 
shock loc.2) 

-29.80 

16.59 

-29.60 

23.60 

3) Turbulent 
(Impinging 
shock loc.2) 

-28.55 

15.34 

-28.26 

21.67 


The magnitudes of the pressure peaks are non-dimensionalized with respect to the 
stagnation pressure of 0.5159 psia. and the magnitudes of heat flux peaks are non- 
dimensionalized with respect to the stagnation heat flux of 20.26 Btu/ft -sec 
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Table D. Coordinates of the triple points and transmitted shock lengths for freestream 
condition set B for the impinging shock location 2 


Case 

1st triple point 
coords, (non-dim) 

2nd triple point 
coords, (non-dim) 

Transmitted shock 
length, mm. 

1) Laminar 

x = -1.70352 

x = -1.08102 

1.03944 

y = -0.24194 

y = -0.54639 

2) Turbulent 

x = -1.77225 
y = -0.25909 

x = -1.11177 
y = -0.57525 

1 .09837 


*** 


All the coordinates are non-dimensionalized with respect to the cylinder radius of 
38.10 mm (1 .5 in) 




Table E. Locations and magnitudes of the pressure and heat flux peaks for freestream 
condition set B 


Case 

Location of the 
pr. peak (deg) 

Magnitude of the 
pr. peak 

Location of the 
h.f. peak (deg) 

Magnitude of the 
h.f. peak 

1) Laminar 
(Impinging 
shock loc.1) 

-34.26 

17.70 

-33.76 

24.63 

2) Laminar 
(Impinging 
shock loc.2) 

-30.14 

15.96 

-29.51 

29.37 

3) Turbulent 
(Impinging 
shock loc.1) 

-32.12 

11.18 

-37.31 

33.15 

4) Turbulent 
(Impinging 

shock loc.2) 

-32.53 

12.67 

-32.45 

49.41 


*The magnitudes of the pressure peaks are non-dimensionalized with respect to the 
stagnation pressure for the undisturbed flow of 7.005 psia, and the magnitudes of the 
heat flux peaks are non-dimensionalized with respect to stagnation heat flux of 
19.78 Btu/ft 2 -sec 
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Fig.2 Problem statement for a hypersonic flow over a cylinder 
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Fig.3 Wall heat flux distribution from LARCNESS computations for a 
laminar flow over a cylinder 




Fig.4 Problem statement for a Type IV shock-shock interference in a hypersonic flow 
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Fig.5 Details of the adapted mesh near the stagnation region of a Type IV shock interference 
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Fig-6 Convergence history of the peak heat flux for the Type IV shock interference 
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Fig.7 Total velocity contours near the stagnation region from the laminar computations for the 
Type IV interference at a freestream Reynolds number of 6,500 
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Fig.8 Pressure contours near the stagnation region from the laminar computations for the 
Type IV interference at a freestream Reynolds number of 6,500 
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Fig-9 First derivative of density contours near the stagnation region from the LARCNESS computations for 
the Type IV interference at a freestream Reynolds number of 6,500 
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Fig. 10 Comparison of the total velocity contours near the stagnation region 
from laminar and turbulent flow computations of Type IV interference 
at a freestream Reynolds number of 6,500 
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Fig.l la Comparison of the wall heat flux distributions from Type IV interference computations at a 
freestream Reynolds number of 6,500 



Fig.1 1b Comparison of the wall pressure distributions from Type IV interference computations 
at a freestream Reynolds number of 6,500 
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Fig.12a Comparison of the wall heat flux distribution from laminar Type IV interference 
computations for two impinging shock locations at freestream Reynolds number of 6 500 



Fig. 12b Comparison of the wall pressure distributions from laminar computations of the Type IV interference 
computations for two impinging shock locations at a freestream Reynolds number of 6,500 
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Fig-14 First derivative of density contours near the stagnation region from the LARCNESS computations 
for the Type IV interference at a freestream Reynolds number of 414,000 
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Fig. 1 5c The total velocity contours near the stagnation region from the turbulent computations of the Type IV 
interference at a freestream Reynolds number of 414,000 
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Fig. 1 6 Streamline contours superimposed on the total velocity contours near the stagnation region from the turbulent 
computation of the Type IV interference at a freestream Reynolds number of 414,000 
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Fig.17a Comparison of laminar and turbulent heat flux distributions from Type IV interference 
computations at a freestream Reynolds number of 414,000 
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Fig.17b Comparison of laminar and turbulent wall pressure distributions from Type IV interference 
computations at a freestream Reynolds number of 414,000 



Fig. 18a Comparison of wall heat flux distributions from Type IV interference computations 
for two impinging shock locations at a freestream Reynolds number of 414,000 
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Fi9-1 8b Comparison of wall pressure distributions from Type IV interference computations 
for two impinging shock locations at a freestream Reynolds number of 41 4,000 
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